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We report the synthesis, crystal structure and magnetic properties of the S=l 2-leg spin- ladder 
compound Na2Ni2(C204)s(H20)2- The magnetic properties were examined by magnetic suscep- 
tibility and pulsed high field magnetization measurements. The magnetic excitations have been 
measured in high field high frequency ESR. Although the Ni(II) ions form structurally a 2-leg lad- 
der, an isolated dimer model consistently describes the observations very well. The analysis of the 
temperature dependent magnetization data leads to a magnetic exchange constant of J = 43 K along 
the rungs of the ladder and an average value of the ^-factor of 2.25. From the ESR measurements, 
we determined the single ion anisotropy to D — 11.5 K. The validity of the isolated dimer model 
is supported by Quantum Monte Carlo calculations, performed for several ratios of interdimer and 
intradimer magnetic exchange and taking into account the experimentally determined single ion 
anisotropy. The results can be understood in terms of the different coordination and superexchange 
angles of the oxalate ligands along the rungs and legs of the 2-leg spin ladder. 



PACS numbers: 75.50.Xx,75.10.Pq, 76.30.-v, 75.30.Gw 

I. INTRODUCTION 

The physical properties of one-dimensional spin sys- 
tems have attracted a lot of attention. For isotropic 
magnetic interactions described in the Heisenberg model 
and low spin quantum numbers the ground state proper- 
ties are strongly influenced by quantum fluctuations and 
often pure quantum ground states are found in a macro- 
scopic system. For the ideal 5= 1/2 antiferromagnetic 
Heisenberg chain (5=1/2 AFHC) with uniform nearest 
neighbor exchange coupling the ground state is a spin 
singlet with gapless excitations 1 ^. Introducing an ad- 
ditional interchain magnetic exchange coupling between 
adjacent spin chains the so-called n-leg spin-ladders are 
formed. These have been discussed very intensively for 
the 5=1/2 case since they exhibit completely different 
ground state properties 3 ^. In particular for even-leg lad- 
ders a non-magnetic spin-liquid ground state with a finite 
energy gap is founcfe. 

For 5=1 already for a single antiferromagnetic chain 
the ground state is a non-magnetic spin-liquid with a 
finite energy gap as proposed by Haldane^ and several 
experimental realizations are known 7 . Only very little is 
known about the ground state and excitations of a S=l 
2-leg spin ladder. From bosonization studies^ and quan- 
tum Monte- Carlo simulations^ a continuous crossover 
from the delocalized Haldane spin gap state to the case 
of localized antiferromagnetic dimers, which also shows 
a spin excitation gap, is predicted. 

Experimentally the situation is complicated by the fact 
that most 5 = 1 systems (usually based on high spin Ni(II) 



transition metal ions) show a strong axial crystal field 
anisotropy D which modifies the ground state of the spin 
systemii. To the knowledge of the authors no realizations 
of an 5=1 spin ladder have been reported so far. 

In this paper we report the synthesis and structure and 
in particular a full experimental and theoretical charac- 
terization of the magnetic properties of a Ni(II) based 
5=1 2-leg spin-ladder material, Na 2 Ni 2 (C204)3(H20)2. 
We performed magnetic susceptibility and high field mag- 
netization studies (up to 55 T) to determine the thermo- 
dynamic properties. In order to study the spectrum of 
spin states in this system, we have used tunable high- 
frequency high-field ESR-spectroscopy at frequencies up 
to 740 GHz in magnetic fields up to 30 T. 

Combining the experimental results from low field 
magnetic susceptibility, high field magnetization and 
ESR we are able to determine all relevant parameters 
of the spin system. In particular we find that the struc- 
turally nearly isotropic 5=1 spin ladder is magnetically 
very well described by isolated spin dimers with an anti- 
ferromagnetic exchange of J = 43 K on the rungs and a 
large axial anisotropy of D = 11.5 K. In order to quan- 
tify the degree of intrachain versus interchain magnetic 
exchange in this system we performed realistic QMC cal- 
culations including the experimentally determined single 
ion anisotropy which support the relevance of the iso- 
lated dimer model for this system. The calculated field 
dependence of the spin states suggests that strong mag- 
netic fields allow to tune the ground state from a non- 
magnetic 5 = to a magnetic 5 = 1 state and at even 
higher fields to the fully polarized 5 = 2 state. This is 
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confirmed by our magnetization data which provide di- 
rect evidence for a magnetic field induced switching of 
the S=0 ground state to an S=l state at B « 30 T and 
to #=2 at B « 60 T . 



II. RESULTS AND DISCUSSION 
A. Synthesis and crystal structure 

The compound NasNis^CUM^O^ (SNOX here- 
after) is synthesized in a hydrothermal reaction from so- 
lutions containing a very high sodium halide concentra- 
tion. Typically Na 2 C 2 4 , Ni(C 2 4 )(H 2 0)2, NaBr and 
distilled water are mixed in a 1:1.9:140:500 molar ratio, 
heated in a sealed autoclave to 225 °C for 5 hours and al- 
lowed to cool slowly to room temperature. The product 
occurs as green crystals of Na 2 Ni 2 (€204)3(1120)2 em- 
bedded in a matrix of solid NaBr. Phase homogeneity 
and purity of SNOX were established by a combination 
of optical microscopy, powder X-Ray diffraction and ele- 
mental analysis. By slowing down the cooling rate, and 
reducing the concentration of Ni 2+ and C2O4 2- we can 
grow single crystals up to 2 mg in lower yielding reac- 
tions. 

Single crystal X-ray structure determination reveals 
SNOX to crystallize in the monoclinic space group P2i/c 
(#14) with a=5.8144(2), b=15.6474(5), c=6.8357(3) A, 
/?=100.56(60) and V=611.37(4) A 3 at 120 K. Refine- 
ment of 118 parameters gave a goodness-of-fit S= 1.074, 
Rl=0.244 and wR2=0.506 on all data, with the largest 
peak (and hole) in the difference Fourier map 0.383 (- 
0.390) eA 3 . 

SNOX is an end member of the series of isostructural 
compounds with the composition Na2T2 (0204)3(1120)2 
with T = Mn(II), Fe(II)A 3 -, Co(II)lA and Ni. Although 
this structure type has been described previously, it 
is useful to recap certain key features here. The Ni(II) 
ion lies on a general position, and experiences a pseudo- 
octahedral coordination environment. It is coordinated 
in a cis geometry by two chelating and crystallographi- 
cally independent oxalate dianions. The remaining coor- 
dination sites are filled by a monodentate oxalate oxy- 
gen atom and a water molecule. Ni-0 bond lengths 
lie within the range 2.020(1) to 2.114(1) A. We note 
that opposite pairs of O have similar Ni-0 lengths, and 
that the 0-0 separations between opposite pairs, for 
which the bond vectors are mostly directed along the 
crystallographic axes, are (d a =4. 163(1), d^=4. 044(1), 
d c =4. 082(1) A along the a, b and c-axes, respectively, 
Fig. 1). The overall structure can be regarded as an 
anionic [N 12(0204) s]^ 1- network with a 1-D ladder like 
topology (Fig. 1), where the metal ions form the vertices 
and the bridging oxalate ions form the bonds. The two 
oxalate ions bridge nickel ions in quite different modes. 
One oxalate forms a symmetric bis-chelating bridging 
mode (with a crystallographic inversion located at the C- 
C centroid) and links pairs of nickel ions (d Ni ... Ni ( rung ^ 



= 5.294(1) A) forming the rungs of the ladder. The sec- 
ond oxalate is unsymmetrically coordinated, chelating to 
one Ni(II) ion and forming a monodentate coordination 
to a second Ni ion, with 1,3-syn anti geometry through 
the bridging carboxylate, this mode forms the ladder legs 
(d-Ni---Ni(ieg) — 5.810(1) A). These structural linkages are 
likely to provide the only significant magnetic superex- 
change pathways. 




Figure 1: Crystal structure of Na 2 Ni 2 (C204)3(H 2 0)2. For 
clarity, Na and H atoms are not shown. The numbers indicate 
the distance (in A) between opposing oxygen atoms on the 
distorted Ni06 octahedra which are tilted by 17 degree with 
respect to the a axis. 



B. Magnetic susceptibility 

We examined several single crystals in magnetic sus- 
ceptibility measurements along the a, b and c axes. The 
measurements were performed using a Quantum Design 
Magnetic Properties Measuring System (MPMS) in ex- 
ternal fields of B ex = 2 T and 5 T in the temperature 
range 2-300 K. In addition we performed measurements 
at different angles between the a axis and the external 
field by rotating the crystal about an axis perpendicular 
to the a axis. 

Fig. 13 shows the temperature dependence of the mag- 
netization along the a, b and c axes in an external field of 
B ex = 2 T. One can clearly see a pronounced maximum 
for all directions with a strong downturn below ~ 50 K 
suggesting a non-magnetic spin-singlet ground state. No- 
tice that the maximum along the a axes (the ladder direc- 
tion) is located at a lower temperature than the maxima 
along the b and c axes (39 K, 47 K and 45 K for the a, 
b and c axes respectively). The difference between the b 
and c axes magnetization is interpreted mainly due to a 
g-factor anisotropy. This behavior in general is expected 
as well for an isotropic two-leg spin ladder as for the two 
limiting cases, an isolated S = 1 Haldane chain or a sys- 
tem of antiferromagnetically coupled dimers^S&^iiS. At 
very low temperatures a small Curie-tail due to a very 
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low concentration of paramagnetic centers (less than 1 
%) masks the approach to zero. 
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Figure 2: Temperature dependence of the magnetization of 
a SNOX single crystal (2.04 mg) at an external field of 
B ex = 2 T along the a (circles), b (diamonds) and c axes 
(triangles) . 

Taking into account the two different exchange path- 
way topologies, we assume a stronger magnetic interac- 
tion J along the rungs than the magnetic interaction K 
along the legs. This leads to an isolated dimer approx- 
imation where each dimer consists of two 5=1 spins on 
the rungs of the ladder. Since all Ni ions on a ladder 
are crystallographically equivalent they share the same 
strength and orientation of the single ion anisotropy D. 
Including the Zeeman energy in the external magnetic 
field, the spin Hamiltonian of the system in the dimer 
approximation is given by 

H = JSiS 2 +^sB(Si + S 2 )+ S ^ S - 

£=1,2 

Here, we assume that the magnetic field is applied along 
the z-axis, the main axis of the crystal field anisotropy 
tensor D_. For pure axial anisotropy this tensor is di- 
agonal, defined as D xx = - 1/3 D, D yy = - 1/3 D and 
D zz = + 2/3 D. To calculate the energy levels for applied 
fields along different directions, we rotate the axis of the 
axial anisotropy tensor by an angle 7 about a perpendic- 
ular axis. Then the spin Hamiltonian can be expressed 
as 

H = JS 1 S 2 +9 7 fJLBB(S lz +S 2 z)+^2 Sig^m^i (2) 

i=l,2 

with the rotation matrix JJ_ and an angle dependent g- 
factor g 7 . 

The isolated dimer approximation may be improved 
by introducing an effective exchange interaction K along 
the legs treated in mean field approximation in the cal- 
culation of the magnetic susceptibility \ using 

Xladder = Xdimer / (1 + KXdimer)- (3) 



To analyse the magnetization data, we developed a 
fit routine which numerical diagonalizes the Hamiltonian 
(0, calculates the free energy A and magnetization M 
for a particular magnetic field orientation 7 and varies 
the parameters in the Hamiltonian to minimize the mean 
square deviation between the data and this model. We 
performed a combined fit of M = x^ex for B ex = 2 T 
along the a, b and c axes using an isotropic magnetic ex- 
change coupling constant J, a fixed absolute value of the 
axial anisotropy of \D\ = 11.5 K (determined by the S=l 
zero field splitting measured by ESR as described below). 
The total susceptibility x consists of the main contribu- 
tion (0 and comprises also a temperature independent 
term xo = Xdia + Xw which includes a diamagnetic 
contribution Xdia and a Van-Vleck paramagnetic suscep- 
tibility Xw of Ni(II), as well as a Curie contribution C/T 
with a Curie constant C owing to paramagnetic impuri- 
ties: 

X = X0 + C/T + Xladder- (4) 

In addition, we used an orientation dependent g-factor g 7 
and, since the orientation of the principal axis (z axis) of 
the magnetic anisotropy tensor is not known, the angle 
7 as parameters for the different orientations. The fit 
results for the B ex = 2 T measurements are shown in 
Fig. [21 as solid lines. 

We obtain an intradimer coupling constant of 
J = 43 K, an interdimer coupling constant of K = K 
and ^-values of g a = 2.206, = 2.291 and g c = 2.256 . 
The angles to the anisotropy axis are j a = 18°, 75 = 90° 
and 7 C = 117°, a Curie constant C = 0.0056 /i^K/dimer 
and a temperature independent contribution xo = 
0.0087 /i#/dimer has been used. The variations of the g- 
values are consistent with the variation of the distances 
to the adjacent oxygens, leading to a larger g- value with 
smaller distance. Using g 7 2 = g|| 2 cos 2 7 + g^ 2 cos 2 7 — 
with g± = gb and g a = 2.206 we find gu= 2.197. The 
knowledge of 7 a , 75 and 7 C allows a determination of the 
orientation of the anisotropy axis in the crystal. 7 a = 18° 
is explained by an anisotropy axis parallel to the line 
connecting both opposite oxygen neighbor ions along the 
a axis (O a axis hereafter). This line is tilted with respect 
to the a axis for all Ni sites (see Fig.QJ right panel) about 
an angle of 17° determined by X-ray crystallography. Re- 
garding 75 we have to consider that with respect to the 
b-axis the O a axis is tilted alternating from ladder to lad- 
der, which leads to an average tilting of 90°, also in very 
good agreement with the fitted value. Finally, 7 C = 117° 
can be explained by taking into account the monoclinic 
structure with Z(a,c) = f3 = 100.6°. The tilting of 117° is 
approximately the sum of the angle between a and c axes 
and an additional tilting of the anisotropy axis with re- 
spect to the a axis towards the c axis of 16°. This is 
10° more than what is actually found for the O a axis. In 
view of the expected accuracy for the determination of 
the anisotropy axis this is reasonable. In summary, from 
an analysis of the angle dependence of the temperature 
dependent magnetization, we can determine the orienta- 
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Figure 3: Relative energies of the spin states, calculated for 
the magnetic field applied perpendicular to the z axis of the 
uniaxial anisotropy tensor: the triplet S = 1 and the quintet 
S — 2 are well separated from the S — ground state by 
an activation energy of roughly |J| = 43 K and 3 |J|, respec- 
tively. Both multiplets exhibit a zero field splitting due to 
crystal field anisotropy. Note the ground state level crossings 
around 30 and 53 T. (for details see text) 



tion of the principal anisotropy axis to be parallel to the 
connecting line of the opposite nearest neighbor oxygen 
ions along the a axis, tilted by 18° with respect to the 
crystallographic a axis mainly towards the b axis. 

The relative energies of the spin states of SNOX calcu- 
lated in the framework of the Hamiltonian @ with the 
parameters yielding the best fit of the magnetization are 
plotted in Fig.|3|for a magnetic field oriented perpendicu- 
lar to the axis of the axial anisotropy tensor. One can see 
that the spin singlet ground state S = is well separated 
from the S = 1 triplet and the S = 2 quintet state. The 
separation energy between singlet and triplet corresponds 
to the coupling constant J, whereas that between singlet 
and quintet corresponds to 3J. The splitting of the ex- 
cited S = 1 and S = 2 states in zero magnetic field is 
induced by the anisotropy D. Both triplet and quintet 
levels split with increasing magnetic field due to the Zee- 
man effect. The role of the zero field and the Zeeman 
splitting for the high field magnetic properties of SNOX 
will be discussed in detail in Section UTTI and Hi Dl 

Note, that analyzing the data with a negative 
D = 11.5 K leads to a similar good agreement between 
the model and the data, resulting in J = 42 K and g- 
values of g a = 2.160, = 2.256 and g c = 2.213 but an 
inconsistent set of angles j a = 90°, 75 = 25° and j c = 46° 
is found. The positive D value is also strongly supported 
by the high field magnetization data described below. 

An appreciable temperature independent contribution 
to the static susceptibility xo can not be ascribed to the 
Van-Vleck paramagnetism Xw alone because typical val- 
ues of Xw fc> r Ni(II) are in the range of 10 -4 emu/molei^ 
leading to 0.0008 /i^/dimer a t B ex = 2T. We assign it 
to a very small concentration of ferromagnetic Ni metal 
clusters, which are formed in the synthesis process as an 
impurity phase. 



C. High Field ESR 

ESR measurements in static magnetic fields up to 15 T 
have been performed with a Millimeterwave Vector Net- 
work Analyzer (MVNA) from AB Millimetre, Parish 
Experiments in pulsed magnetic fields up to 30 T have 
been carried out at the National Laboratory of Pulsed 
Magnetic Fields, LNCMP, Toulouse, France, using a 
Fabry-Perot cavity optically pumped by a CO2 laser for 
generation of microwaves. Details of both experimen- 
tal set-ups can be found in Ref. For the ESR mea- 
surements a powder sample of SNOX was prepared as a 
pressed pellet with a small amount of a diphenyl-picryl- 
hydrazyl (DPPH) marker attached for calibration of the 
magnetic field strength. 

At a base temperature of 4.2 K no ESR signals are 
found which can be associated with the resonance of bulk 
Ni(II) spins in the sample. However, it is possible to 
assign four modes to the bulk ESR spectrum of SNOX 
which develop in the spectrum with increasing temper- 
ature (Fig. 21) • These modes are numbered from 1 to 4 
in order of their appearance in ascending magnetic field 
and marked by arrows in the inset of Fig. 21 The position 
of the modes in a magnetic field B\ es depends on the ex- 
citation frequency v as would be expected for magnetic 
resonance. Respective resonance branches v versus B\ es 
are plotted in the main panel of Fig. 21 In the ESR exper- 
iment in static magnetic fields with v < 550 GHz mode 1 
is the strongest in intensity whereas modes 2, 3 and 4 ap- 
pear as weak satellites on the nonlinear background. At 
a fixed excitation frequency v = 300 GHz the intensity 
/ of mode 1 has been studied as a function of temper- 
ature. The result is shown in Fig. [51 I(T) exhibits an 
activation temperature dependence similar to that of the 
static magnetization M(T) in Fig. [3 Because the ESR 
intensity is determined by the static susceptibility x °f 
the spins participating in the resonance 12 the common 
temperature behavior of I{T) and M(T) = x(T)B gives 
evidence that in the ESR experiment the same spins are 
probed which contribute to the static magnetic proper- 
ties. Mode 1 has a maximum intensity at T ~ 30 — 35 K 
which roughly corresponds to the energy separation be- 
tween the singlet ground state and the first S = 1 ex- 
cited state (Fig.[3J). This allows to unambiguously assign 
this mode to a resonance excitation within the S = 1 
multiplet. Owing to the weakness of modes 2, 3 and 4 
they can be detected only in a limited temperature range 
~ 30 — 60 K. At smaller temperatures they are not visi- 
ble due to the strong decrease of their amplitude whereas 
at higher temperatures they also broaden. Nevertheless 
their occurrence in the temperature window where the 
static magnetization is strongest (Fig. |2| allows to assign 
them to the bulk ESR response of the SNOX sample too. 

A remarkable feature of the frequency /field diagram 
in Fig. 21 is the intercept of the strongest ESR branch 1 
with the frequency axis at v§ « 239 GHz. It implies a 
finite energy gap A = v^h/ks = 11.5 K for this reso- 
nance excitation. This zero field gap can be straightfor- 
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Figure 4: Main panel: Relationship z/ versus B\ es for the 
ESR modes i = 1 ... 4. Solid line is a theoretical fit of the 
data. Dashed lines is a linear approximation of branches 2, 3 
and 4 to high fields. Gray error bar indicates a field region 
where a broad absorption mode is observed in the pulsed field 
ESR measurement. Inset: ESR spectrum in a static (bottom 
curve) and pulsed magnetic field ESR experiment (top curve). 
Sharp lines are due to the DPPH field marker. Resonance 
modes associated with a bulk ESR spectrum of Ni-oxalate are 
indicated by black arrows and numbered. The anticipated 
positions of modes 2, 3 and 4 at high magnetic fields are 
marked by gray dashed arrows, (for details see text). 
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field B ex . As a result orientations perpendicular to the 
z axis of the uniaxial anisotropy tensor will dominate in 
the spectrum whereas other directions may form a broad 
nonlinear background. Thus for a semiqualitative under- 
standing of the field dependence of the resonance modes 
it is sufficient to consider the energy diagram for B ex JL z 
shown in Fig. 03 The field dependence of branch 1 can 
be very well reproduced (solid line in Fig. |3} assuming 
that the resonance transition occurs between the | — ) 
and | + ) states of the excited triplet (Fig.EJ). In fact, this 
usually " forbidden" transition^ has the strongest inten- 
sity in the spectrum for v < 550 GHz. This is due to 
the strong anisotropy of the crystal field yielding a large 
zero field gap A and mixing strongly the pure spin states 
|— 1), |0) and |+1) into the | — ) and | + ) states for 
B ex JL z in the field regime g\iB^ex < A. Indeed, in a 
pulsed field ESR experiment the relative intensity of this 
mode becomes appreciably smaller (see Fig.0] inset). For 
the excitation frequency of 739 GHz the resonance field 
amounts to B\ es = 11.2 T. For this field strength not the 
crystal field but the magnetic field direction determines 
a quantization axis. Thus the mixed | — ) and | + ) states 
transform into almost pure | — 1 ) and | +1 ) states which 
do not fulfill the ESR selection rule AS Z = ±1. 

In contrast to mode 1, the exact assignment of modes 
2, 3 and 4 is difficult owing to their small intensity. The 
relationship v versus B[ es has a steep slope because both 
energy levels, | — ) and | + ), depend on magnetic field di- 
verging from each other with increasing B ex . The much 
smaller slope of branches 2, 3 and 4 suggests that reso- 
nance transitions within the neighboring pairs of states, 
like e.g. (|-) <-> |0)), (|0) <-> | + )) etc., of the multi- 
plets 5 = 1 and S = 2 should be involved. Branch 2 has 
a much smaller frequency offset as compared to A and 

can be associated with the transition | ) <-> | — ) in 

the quintet 5 = 2. Branches 3 and 4 might be gapless, 
which is the case for transitions | — ) <-> | ) in triplet 
as well as in quintet. Unfortunately, it was not possible 
to achieve a better resolution of modes 2, 3, and 4 in 
the pulsed field measurements. Only a broad absorption 
mode is observed in the field region where resonances 2, 
3 and 4 are anticipated (Fig. EJ). 



Figure 5: Intensity of the strongest ESR mode 1 as a function 
of temperature. Note an activating behavior similar to the 
static magnetization M(T) (Fig.0. 



wardly identified with the zero field splitting of the first 
excited triplet in Fig. 03 Numerical solution of the Hamil- 
tonian © yields A = \D\. Therefore, from the ESR data 
one obtains an accurate estimate of the absolute value of 
the single ion anisotropy parameter \D\ = 11.5 K which 
is essential for the understanding of the static magnetic 
properties of SNOX as discussed in Section III Bl and Til Dl 
ESR measurements of a powder sample imply the av- 
eraging of all possible orientations of the microcrystal- 
lites in a powder with respect to the applied magnetic 



D. High field magnetization 

The high field magnetization measurements were car- 
ried out at the LNCMP, Toulouse, as well. Measure- 
ments of the static magnetization have been performed 
in pulsed magnetic fields up to 55 T with a pulse dura- 
tion of 250 ms (raising time 40 ms) by an inductive tech- 
nique. A system of two concentric pick-up coils in oppo- 
sition to each other equipped with additional compensa- 
tion coils was utilized. The measured signal dM/dB was 
integrated numerically. The sample temperature can be 
varied between 1.4 K and 300 K. 

As can be seen from the energy diagram of the spin 
states of SNOX (Fig. the splitting of the states in a 
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magnetic field yields a level crossing of the ground state 
with the lowest triplet state | — 1 ) at a field Bci and a 
second level crossing of the | — 1 ) state with the | — 2 ) 
quintet state at field Bc2 with Bci < Bc2- This leads 
to a step-like behavior in high field magnetization mea- 
surements. If the field is applied along the z axis, the 
critical field Bci can be estimated using the equation 
Bci = (J — D/3)/(gfiB/kB), depending strongly on 
the sign of D. For a field perpendicular to the local 
anisotropy axis, a negative D (leading to higher criti- 
cal fields for B || z) pushes the critical field to lower 
fields and vice versa for a positive D. Since in powder 
measurements the perpendicular situation dominates the 
spectrum of the spin states, the field dependent magne- 
tization, among others, can determine the sign of the 
anisotropy D. 

The critical field strengths Ba depend on the magnetic 
exchange J and an exchange energy of 43 K corresponds 
to roughly 30 T for Bci- Therefore high magnetic fields 
are needed to observe the magnetization steps. We per- 
formed measurements at several temperatures (1.47 K, 
4.2 K, 10 K and 26 K) on a powder sample with a mass 
of 23.2 mg in magnetic fields up to 55 Tesla. The results 
are shown in FigEl In the low temperature experiments 
the first magnetization step is clearly seen, whereas it 
smoothes at higher temperatures due to the thermal av- 
eraging process. The critical field Bci can be obtained 
from the derivation of the magnetization dM/dB which is 
shown in FigQfor T = 1.47 K and T = 4.2 K (inset). At 
both temperatures a sharp peak is observed at 29 Tesla. 
Additionally a second increase of the magnetization at 
much higher fields of about 50 Tesla can be anticipated. 
This step is also observed strongly broadened in the 10 K 
measurement. At 26 K a linear field dependence of the 
magnetization is found. 

The solid lines in FigEl and FigQ describe simula- 
tions using the dimer model with the calculated g val- 
ues g\\ (see above) and g± = ^ in a powder average 
M = (M||+2 Mj_)/3 . To describe the high field mag- 
netization well we need to use a coupling constant of 
J = 44 K instead of J = 43 K. In particular the dif- 
ferential magnetization is very well reproduced by the 
simulation, describing as well the width as the center po- 
sition of the peak. A small deviation is found in the field 
range of 30 - 33 T in the T = 1.47 K measurement, which 
may depend on a little shift of the sample induced by the 
external field. 

The simulations for a negative D = -11.5 K did not 
describe the data reasonably well for any J. In this case 
the first magnetization step is consistent with J = 43 K, 
whereas the second step is predicted at too low a field. To 
reproduce the second step well one has to use J = 45 K 
which results in a too large a value of Bci- 

From Fig0 we can determine the critical field Bci 
of the powder sample to Bci = 29 T. In our simula- 
tions, we get Bci = 27.7 T for B parallel to the mag- 
netic anisotropy axis and Bci = 30.7 T for B perpen- 
dicular to the magnetic anisotropy axis, resulting also 




Magnetic Field (T) 

Figure 6: High field magnetization of a SNOX powder sample. 
The figure shows measurements at 1.47 K and 4.2 K, the inset 
measurements at 10 K and 26 K. The solid line represents the 
simulation as described in the text. 




-10 10 20 30 40 

Magnetic Field (T) 



Figure 7: Differential magnetization curves dM/dB for tem- 
peratures T = 1.47 K and T = 4.2 K (inset). The open circles 
are calculated from measured data, the solid lines represent 
simulations as described in the text. 



in Bcio — 29 T in a powder average which is shown in 
the figure. Notice, that the powder averaging leads to 
a broadening of the magnetization steps. For the second 
step, we get Bc2 = 67.1 T and Bc2 = 51.3 T for B paral- 
lel and perpendicular, respectively. This difference leads 
in a powder sample (average over all angles, which we 
calculated for angles from 0° to 90° in steps of 1° with 
varying g from gu to g±) at T = 1.47 K to a continu- 
ous double-sigmoid like increase of the magnetization in 
the field range of 50 - 70 Tesla, of which only the first 
increase is observed in the measurement. 



E. Quantum Monte Carlo Calculations 

To perform a quantitative analysis of the ratio 
R=J/(K+|J|) of the intra- to inter-dimer coupling, we 
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have compared the observed magnetic susceptibility and 
the magnetization with Quantum Monte Carlo (QMC) 
calculations. These calculations explicitely include the 
effects of the finite single ion anisotropy D = 11.5 K. 
This is different from previous QMC studies of spin-1 
ladders 10 . The QMC method we employ is the stochas- 
tic series expansion (SSE)A2iSft, which directly samples an 
expansion of the partition function in powers of 1/T until 
convergence. Within the SSE the cumulants of this power 
series are evaluated using diagonal and non-diagonal up- 
dates within the space of exchange operators from eqn. 
(0). The latter update is greatly improved by introduc- 
ing directed loop- updates^. We have carried out SSE 
calculations on systems with up to 200 x 2 sites and in- 
verse temperatures (3 < 10. The system sizes have been 
checked to represent the thermodynamic limit within the 
temperature range studied. In all figures depicted the 
statistical QMC errors are much smaller than the sym- 
bol sizes and have been omitted. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
Reduced Temperature T/(J+K) 



Figure 8: Quantum Monte Carlo simulations for R — 1 (open 
circles), 0.95 (open triangles) and 0.9 (open squares). We 
plotted the reduced susceptibility versus reduced tempera- 
ture. The solid circles show the measured data for B — 2 T 
along the a axis corrected by a Curie tail and a temperature 
independent contribution. The inset shows the high temper- 
ature results. 

Our calculations show a significant change of the 
shape of the susceptibility at low temperatures as a 
function of R. In the plot of the reduced suscepti- 
bility Xred = X (J+K) versus reduced temperature 
T re d = T/(J+K) the peak increases and narrows with 
decreasing R down to R = 0.7 and then decreases and 
broadens with further decreasing The calculations 

have been performed for D/J = 0.27, which corresponds 
to D = 11.5 K for J = 42.6 K, similar to the fit 
results described above. The results for R = 1, 0.95 
and 0.9 are shown in Fig|S| compared with the scaled 
experimental data. For this scaling we converted the 
magnetization data for B = 2 T parallel to the a axis 
to the reduced representation and vary the ^-factor and 
the sum of coupling constants J+K (which scales T re d 



and Xred) to fit the data to the high temperature range 
(2< T re d<8) of the QMC calculations. Notice, that the 
scaled experimental dataset appreciably depends on the 
^-factor (in the range of 2.2 to 2.25). Nevertheless, all 
possible transformations lead to a value of R w 0.95 - 
1.00 for total magnetic exchange strenghts between 43 K 
(R = 1) and 45 K (R = 0.95). The solid circles in FigEI 
show the best result of the transformation using a total 
coupling constant J+K = 45 K, a ^-factor g = 2.23, a 
temperature independent term of 0.0085 /i# /dimer and 
a Curie constant of 0.0045 /i^K/dimer. A comparison 
with the QMC data in the low temperature regime 
of T re d<1.5 gives an R of 0.95, which results in an 
intradimer coupling of J = 42.75 K and an interdimer 
coupling of K = 2.25 K. 

QMC calculations for the high field magnetization of a 
weakly coupled dimer system were done for a reduced 
temperature T/J = 0.034 and a reduced anisotropy 
D/J = 0.27 corresponding to a sample temperature 
T = 1.47 K and D = 11.5 K if J + K = 43 K. The 
calculations for R = 1, 0.95 and 0.9 are compared with 
the measured data in FigO The measured data were 
scaled by setting the magnetization at B ex = 40 T to 
M = 0.5 and using a value J + K w 43 K. A variation 
of J + K mainly shifts the position of the step and does 
not affect the broadening of the step. The QMC calcu- 
lations show a broadening of the magnetization step for 
decreasing R, i.e. increasing magnetic exchange strength 
along the legs of the ladder. The QMC calculations are 
possible only for B parallel to the anisotropy axis up to 
now. Therefore, comparing with the powder measure- 
ments one has to consider that the magnetization step in 
the experimental data set is additionally broadened by 
the presence of the axial anisotropy D. Therefore, the 
step gradient of the appropriate QMC calculation has to 
be sharper than the measurement performed on a pow- 
der sample. For this reason, the experimental data are in 
agreement only with the calculation for R — 1, i.e. the 
uncoupled dimer model. 



III. THE ANISOTROPY OF THE MAGNETIC 
EXCHANGE STRENGTH 

An unexpected result of the present study is the very 
large difference between the strength of the magnetic ex- 
change interaction along the rungs (J = 43 K) and legs 
(\K\ < 2 K) of the spin ladder. This can be qualita- 
tively understood considering the different coordination 
and superexchange angles of the exchange mediating ox- 
alate molecules. 

On the rungs the oxalate forms a \i- 1,2, 3, 4 bridge be- 
tween two Ni ions. This bridge provides two symmetric 
superexchange pathways. On each path the Ni 3d x 2 -y 2 
orbitals have an enhanced electron probability density ex- 
tending directly towards the corresponding oxygen ions. 
There is a direct overlap of the Ni 3d x 2 -y 2 and the O 
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Reduced Field B/(J+K) 

Figure 9: QMC calculations of the field dependent magneti- 
zation for different ratios of R with rescpected D/J — 0.27. 
Open circles: R — 1, open triangles: R — 0.95, open 
squares :R = 0.9. The calculations were done for a reduced 
temperature T / J — 0.034, corresponding to J + K — 43 K 
and T — 1.47 K. The solid squares are the normalized mea- 
surements scaled with J + K = 43K 




Figure 10: Ni-(ox)-Ni coordination and proposed magnetic 
exchange pathways on the rungs (left panel) and on the legs of 
the spin ladder (right panel). The involved Ni 3d and oxygen 
2p orbit als are shown in black. 



O 2p orbit als involved in the Ni-0 bonds are not over- 
lapping with each other and therefore either the carbon 
atoms or orthogonal O 2p orbitals are involved which 
strongly suppresses the strength of the superexchange 
mechanism and may even lead to a weak ferromagnetic 
exchange. Similar /i-1,2,3 oxalato bridges are reported 
for Cu(II) systems with magnetic exchange strengths in 
the range of -0.2 to 0.3 K2i2£ consistent with the esti- 
mate of \K\ < 2 K for SNOX. 



IV. CONCLUSION 

In this paper we report the synthesis, crystal structure 
and magnetic properties of the structurally well defined 
5 = 1 spin ladder compound Na2Ni2 (€204)3(112 0)2. We 
performed magnetic susceptibility and pulsed high field 
magnetization measurements to examine the thermody- 
namic properties of the spin system. The magnetic exci- 
tations have been examined in high field high frequency 
ESR. Although the Ni(II) ions form a structural ladder, 
an isolated dimer model on the rungs describes the ob- 
servations very well. The analysis of the temperature 
dependent magnetization data leads to a magnetic ex- 
change constant of J = 43 K along the rungs of the lad- 
der. We deduce the principal axis of the axial anisotropy 
tensor. From ESR measurements, we determined the sin- 
gle ion anisotropy to D = 11.5 K. To prove the validity 
of the isolated dimer model, we compare the single crys- 
tal magnetic susceptibility measurements with Quantum 
Monte Carlo calculations for a S=l spin ladder, which 
were carried out for several ratios between interdimer and 
intradimer exchange and take into account the experi- 
mentally determined single ion anisotropy D. It results 
in a upper limit of 5 % interdimer exchange coupling. 
This strong anisotropy of the magnetic exchange inter- 
action can be traced back to the different superexchange 
pathways on the oxalate molecules along the rungs and 
legs of the spin ladder. 



2p wave functions forming a bonds. The polarized 2p 
orbitals themselves are strongly overlapping. Therefore 
the intermediate carbon atom is not involved in the su- 
perexchange mechanism resulting in a strong antiferro- 
magnetic superexchange interaction (see fig. HU|) . In the 
literature several /i-1,2,3, 4 oxalato-bridged Ni(II) dimer 
and chain systems are reported (see e.gffi and refer- 
ences therein) with J values in the range of 20 to 42 K. 
The value of 43 K in SNOX is at the upper limit of the 
reported range. 

For the oxalate bridge along the legs of the spin ladder 
the situation is completely different. In this case a /i-1,2,3 
oxalato bridge is formed. As shown in the right panel of 
fig. El on ly one Ni ion forms two covalent bonds (one 
with the x 2 -y 2 orbital and one with the 3z 2 - r 2 orbital) 
with the oxalate molecule whereas of the second Ni ion 
only the 3z 2 - r 2 orbital forms one covalent bond. The 
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